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A numerically implementable Multi-scale Many-Body approach to strongly correlated electron 
systems is introduced. An extension to quantum cluster methods, it approximates correlations 
on any given length-scale commensurate with the strength of the correlations on the respective 
scale. Short length-scales are treated explicitly, long ones are addressed at a dynamical mean-field 
level and intermediate length-regime correlations are assumed to be weak and are approximated 
diagrammatically. To illustrate and test this method, we apply it to the one dimensional Hubbard 
model. The resulting multi-scale self-energy provides a very good quantitative agreement with 
substantially more numerically expensive, explicit Quantum Monte-Carlo calculations. 



I. INTRODUCTION 

Strongly correlated electron systems are characterized 
by a manifold of complex, competing phenomena, which 
emerge in the thermodynamic limit. The underlying 
mechanisms involve correlations on all length scales. 
Currently no single feasible numerical method exists to 
accurately address these correlations at all lengths. Both 
finite size and mean-field calculations alike are faced with 
these limitations. However, interest in areas such as 
quantum phase transitions and magnetically-driven su- 
perconductivity point to the need for numerical schemes 
that accurately bridge the short and long length scales. 

Quantum cluster techniques^ , constitute a good staring 
point in addressing the entire range of correlations by 
dividing the problem into two length regimes; explicitly 
solving for short ranged correlations and approximately 
for the remaining longer length-scales. One such method, 
the dynamical cluster approximation (DCA)^, maps the 
lattice problem onto an embedded cluster problem. In 
doing so, short-ranged correlations within a given cluster 
are treated accurately, while the remaining longer-ranged 
correlations are approximated on a dynamical mean-field 
level. However, the extent of length-scales which can 
thus be accurately addressed beyond the mean-field level 
is severely limited by the numerical expense involved. 

This limitation results in an inadequate treatment 
of medium ranged correlations which are outside the 
scope of explicit calculations. We introduce a Multi- 
scale Many-Body (MSMB) approach which addresses 
each length-scale using approximations adequate for the 
strength of the correlations on the respective scale. The 
strongest, local and short ranged correlations are well ac- 
counted for in traditional, numerically exact implementa- 
tions of the DCA. Correlations, except for in the vicinity 
of phase transitions, fall off rapidly with distance and 
are hence considered weaker in the intermediate length- 
regime. However, these correlations remain significant 
and will hence be approximated diagrammatically. Only 
the remaining third regime of the longest length scale will 



be treated at the dynamical mean-field level. 

The perturbative inclusion of correlations on an inter- 
mediate length-scale within a multi-scale approach has 
previously been explored by Hague et. aA Contribu- 
tions to the single-particle self-energy on various scales 
were linked in a hybrid approach. However, the inher- 
ent perturbative nature of the approach limited it to 
high temperatures and/or weak coupling strengths. In 
this work, we present a non-perturbative two-particle di- 
agrammatic approach to the intermediate length-regime. 

To illustrate and test this scheme, this MSMB is ap- 
plied to the one dimensional Hubbard model"*. While a 
formally simple model, the Hubbard model contains a 
multitude of the underlying physics of correlated elec- 
tron systems. It therefore lends itself ideally as a bench- 
mark for the method. We show that the MSMB approach 
yields results in very good quantitative agreement with 
explicit large cluster calculations. 

Before we proceed to provide detailed results in sec- 
tion we first establish the theoretical basis for the 
method. A subsequent outlook on further developments 
in MSMB techniques is provided in section IVTl conclud- 
ing with a brief summary. 

II. FORMALISM 

For simplicity, we will use the one dimensional Hub- 
bard model to illustrate the MSMB formalism. This low 
dimension is also the most difficult regime for quantum 
cluster approaches like the DCA. The Hubbard Hamilto- 
nian is given by 

H = - t{c\^Cj„ + H.c.) + e ^ rij^ 

<ij> id 

+C/^(n,|-l/2)(na-l/2) 

i 

with creating an electron of spin a at site i and local 
density riic = cl^Cia. The first part, the kinetic term, al- 
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lows hopping between adjacent lattice sites with transfer 
integral t. The second term is the on-site Coulomb repul- 
sion making a doubly occupied lattice site unfavorable. 
Throughout the remainder of this paper we choose the 
bare bandwidth W — At as the unit of energy by setting 
t — 0.25 and work at fixed filling n — 0.75. 

A. DCA 

The DCA is a systematic quantum cluster theory that 
maps the lattice problem onto a self-consistently embed- 
ded cluster problem. It is an extension of the dynam- 
ical mean-field theory (DMFT)^-® which systematically 
incorporates non-local correlations. In the limit when 
the cluster size is one (i.e. single site), it recovers the 
purely local DMFT solution, systematically incorporates 
non-local corrections as the cluster size increases, and fi- 
nally becomes exact when the cluster size equals the size 
of the lattice. 

The respective approximations for the DMFT and 
DCA may be derived by approximating the Laue function 
which describes momentum conservation at the vertices 
of the irreducible diagrams: 

A, , — \^ „ir(ki+k2-k3-k4) /on 



^exact — N^ki+k2M + ki ■ (3) 

In the DMFT, the Laue function is approximated with 
^DAiFT — 1 for all combinations of fci, k2, and k^. 
In doing so, all electron propagators in the self-energy 
diagrams may be averaged over the first Brillouin Zone 
(BZ); thus relinquishing any momentum dependence of 
the self-energy. Hence, the DMFT lattice Green's func- 
tion contains local correlations of the system but is un- 
able to capture non-local correlations. The DCA sets out 
to systematically include these non-local contributions. 
This is accomplished by partially restoring momentum 
conservation of the irreducible vertices. We divide the 
BZ into Nc identical discrete sub-cells as illustrated in 
Fig. n The center of each cell is labeled by if, and the 
surrounding points by fc, so that any arbitrary k — K+k. 
In the DCA, this partial momentum conservation is ex- 
pressed by the Laue function: 

A_DCA — Nc5ki+K2.K3+K4 ■ (4) 

Therefore, all propagators may freely be summed over 
intra-cell momenta k, yielding the coarse grained Green's 
function 

G{K,iLOn) = ^ ^ G{K + k,iuJn) ■ (5) 
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FIG. 1: In the DMFA all electron propagators are averaged 
over the entire Brillouin Zone (BZ); effectively mapping the 
lattice onto a single point (top). In the DCA, we break the 
BZ into several sub-cells which are now in turn averaged over 
mapping the lattice onto a finite sized cluster (bottom). 




FIG. 2: Self-consistency loop for the DCA. 



In so doing, only momentum conservation of magnitude 
Afc < {2-k/Nc) is neglected, while larger inter-cell trans- 
fers are preserved. The resulting self-energy diagrams are 
now those of a finite cluster of size where each lat- 
tice propagator has been replaced by its coarse-grained 
analog, and the remaining cluster problem is defined by 
G{K,iujn)- We can write for the DCA lattice Green's 
function 



iUn + ti~e{K + k)-^{MiK + k),iUn) 

(6) 

where M{k) is a function which maps momentum k re- 
siding in a certain sub-cell of the BZ to its cluster mo- 
mentum K and the lattice self-energy is approximated 
by that of the cluster problem. 

The remaining embedded cluster problem must be 
solved with a self-consistency requirement that the 
Green's function calculated on the cluster Gc{K,iu!n) = 
G{K,iuJn)- Fig. 13 depicts the corresponding DCA algo- 
rithm: Starting with an initial guess for the self-energy, 
we construct the coarse-grained Green's function G from 
the corresponding lattice G (see Eq. EJ. In the next 
step, we utilize one of the many available cluster solvers 
to determine the cluster self-energy. This is the nu- 
merically most involved step and a variety of numeri- 
cal techniques may be applied. At this point, we use 
the new estimate for the self-energy to re-initialize the 
self-consistency loop. It is important to notice that in 
this procedure only the irreducible lattice quantities are 
approximated by their cluster equivalent - i.e. the self- 
energy. 

The DCA has been successfully implemented with a va- 
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riety of cluster solvers of which some are exact but limited 
in cluster size, while others are applicable up to larger 
length-scales but involve varying degrees of approxima- 
tion. Some of the cluster solvers which have been used 
in conjunction with the DCA include the non-crossing 
approximation (NCA)i, the fluctuation exchange ap- 
proximation (FLEX)i2 and the Quantum Monte Carlo 
(QMC)2 method. While NCA and FLEX involve vari- 
ous levels of approximations, QMC is of special interest 
since it provides an essentially numerically exact solution 
to the problem. Although the QMC constitutes a precise 
cluster solver, it becomes prohibitively expensive for large 
clusters. The range of applicability of exact calculations 
is thus restricted to relatively short length-scales. How- 
ever, various properties of strongly-correlated systems 
are not accounted for (e.g. Mermin- Wagner theorem^) 
due to the absence of long-ranged fluctuations in these 
solutions. 



B. Multi Scale Method 

The inability of a single solver within the DCA to 
numerically address long-ranged correlations explicitly, 
motivates a MSMB approach where the problem is di- 
vided further, incorporating a third, intermediate length- 
regime. In this approach, the lattice problem is mapped 
onto two clusters of different size each of which con- 
tributes correlations of length-scales up to the linear ex- 
tent of their respective cluster size. The respective cluster 
problems are addressed using approximations adequate 
with the strength of the correlations on the respective 
scale. We choose a small DCA cluster of size nP to 
be solved using the QMC, thus explicitly accounting for 
the shortest ranged correlations in the system. Next, we 
invoke a second, larger cluster of size Nc to address 
the intermediate length regime. Except in the vicinity 
of phase transitions, correlations on these longer length- 
scales are weaker and the corresponding self-energy is 
approximated diagrammatically. 

We build a suitable approximation by considering the 
single-particle self-energy (the Hartree term is not explic- 
itly shown) written in terms of the reducible vertex F in 
both the particle-hole and particle-particle channel as de- 
picted in Fig. 13 _F in turn is related to the irreducible 
vertex F via the Bethe-Salpeter equation (see Fig. ^ 

F{k, k' , q; lUJri, i^^'n, iv) = T{k, k' , q; iw„, zw^, iv) 

+ T{k, fc", q; luj^, iJ^, iv) x yi'ik" , q; icj„, tv) 
X F{k",k',q;iujl[,iuj'„,iv) 

(7) 

In perturbation theory, the approximations to the self- 
energy are often made at the level of the irreducible ver- 
tex function. In the simple approximation F = [/, the 
resulting self-energy diagrams in Fig. |3| are those of the 
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FIG. 3: Diagrams relating the self-energy to the reducible 
two-particle longitudinal spin and charge (top) and transverse 
spin {middle) vertices. A similar relation is obtained for the 
particle-particle channel {bottom). 
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FIG. 4: Bethe-Salpeter equation relating the reducible two- 
particle vertex F to the irreducible vertex function F. 



FLEX. The resulting method provides qualitatively cor- 
rect long-range properties, but short ranged correlations 
are addressed inadequately - e.g. local moment forma- 
tion. 

The failure of the FLEX, and similar perturbative ap- 
proaches based on the irreducible vertex, occurs in the 
stronger coupling regime. This is a result of neglecting 
higher order correction terms which re-normalize the ver- 
tex. In a non-perturbative approach, Y.M. Vilk et.al. in 
Ref.T^ considered a renormalized static vertex to address 
this problem. 

In the MSMB method however, we introduce the fol- 
lowing non-perturbative approximation for the interme- 
diate length-regime: The large cluster irreducible two- 
particle vertex is approximated by the small cluster ir- 
reducible vertex thus preserving the exact short ranged 
correlation (on the small cluster) and approximating the 
intermediate ranged ones (on the large cluster): 

(8) 
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where momenta labeled with subscript 1 are cluster mo- 
menta on the small cluster, while subscript 2 denotes 
large cluster momenta. The approximated self-energy 
on the large cluster is evaluated from the Dyson equa- 
tion depicted in Fig. |2| It preserves all short ranged 
correlations and includes long-wave length contributions 
which emerge from the inclusion of large cluster correc- 
tions through the bare bubble x'^, i-e. the reducible part 
of the vertex. However, by implementing the DCA on 
a two-particle irreducible level (i.e. replacing the large 
cluster irreducible vertex with the small cluster equiva- 
lent), the resulting self-energy is inherently different from 
the single-particle DCA self-energy that one would obtain 
from the usual, direct DCA/QMC simulation of the large 
cluster (see Sec. Ill A|l . Only in the case of infinite dimen- 
sions, where the problem becomes purely local, are these 
two approaches equivalent. 

A few remarks about the above approximation are in 
order. The approximation in Eq. |S1 breaks the crossing 
symmetry of F (as it is related to F by Eq. [7|), since F 
now contains long-ranged corrections beyond the linear 
extent of the small cluster only in the channel it is calcu- 
lated in. Hence, the self-energy contribution have to be 
evaluated in all three depicted channels (see Fig.|31). Fur- 
thermore, this MSMB approach requires a full knowledge 
of the momentum and frequency exchange on the cluster 
and hence the evaluation of F on the small cluster still in- 
volves extensive numerical calculations and is limited by 
storage/memory requirements (see Sec. I VII for details). 

This difficulty of obtaining the small cluster irreducible 
vertex F, necessitates a further long-ranged approxima- 
tion for F. In ref.ll4l Abrikosov et. al. have, in their study 
of Fermi liquid theory, identified a sub-class of self-energy 
diagrams which convey the long length-scale properties of 
a system. These are constructed from irreducible vertex 
functions with zero external momentum and frequency 
transfer. In applying this restriction we expect to cap- 
ture long ranged correlations which are characterized by 
small momentum transfers q. The resulting long wave- 
length approximated irreducible vertex, denoted by the 
index A, is given by 

T{Ki,K[,Qi]iujn,iuj'^,ii'ri) ~^ 

F^(i4:i,/f(,(3i = 0;icj„,icj^,w„ = 0) . 

(9) 

The MSMB approximation at the level of the irre- 
ducible vertex is expected to provide the best description 
for the most dominant two-particle processes. Further- 
more, any phase transitions would manifest themselves 
in instabilities of the reducible vertex which contains the 
singular structure while the irreducible vertex would re- 
main analytic throughout. However, the MSMB result is 
expected to remain most valid away from phase transi- 
tions where intermediately ranged correlations are weak. 
When these correlations are not weak, a diagrammatic 
approximation of the self energy may contain significant 
errors, but such an approach can nonetheless represent 



non-trivial aspects of strong correlations such as non- 
Fermi-liquid behaviortiMS,. Further details on the calcu- 
lation of the A self-energy are provided in section Hill 

The self-energy obtained in this A-approximation only 
correctly accounts for long-ranged fluctuations of the sys- 
tem. The remaining short length-scale contributions have 
been neglected and have to be accounted for separately 
as will be discussed in a subsequent section 111171 



C. A Conserving Approximation 

There are currently two widely used approaches which 
are sufficent to show an appr oximat ion is conserving. For 
one, G. Baym (see Ref. has shown that, with lo- 
cal conservation of spin, charge, momentum and energy 
at each vertex, it is sufficient for the irreducible vertex 
to be a functional derivative of the large cluster self- 
energy T{G{K,uj),U) = 5'L{G{k,uj),U)/5G{k,uj). How- 
ever, within both the DMFT and DCA momentum con- 
servation is partially violated at each of the internal 
vertices as described by their respective Lane functions 
(Eq. I^J. This violates certain Ward identities and hence 
neither the DMFT nor the DCA constitutes a conserving 
approximation (see ref. for details). 

In an alternative approach, following the arguments 
of Baym and Kadanoff^'^ an approximation can also be 
shown to be conserving as long as it fulfills both the re- 
quirements of 1) the inversion symmetries of F to be 
preserved and 2) that the two-particle correlation func- 
tion tends to the single particle Green's functions via the 
Dyson's equation. It is straightforward to show that the 
approximation for the self-energy on the large cluster sat- 
isfies these requirements. Thus, the MSMB formalism 
based on the explicit small cluster vertex (Eq. |S| only) 
constitutes a conserving approximation for the large clus- 
ter as it restores momentum dependence of the large clus- 
ter self-energy. 

With the introduction of the A-approximation (Eq. O 
short ranged correlations are neglected and have to be 
supplemented as discussed in the subsequent section. 
The resulting self-energy breaks the conservation laws 
with corrections of order 1 /Lc (where Lc is the linear 
size of the larger cluster). It should be reiterated how- 
ever that the necessity of the second (long wavelength) 
approximation is only temporary until the hurdle of large 
memory requirements can be met. 



D. Ansatz 

To account for the omitted set of short length-scale self- 
energy diagrams in the A-approximation, we substitute 
appropriate diagrams from the small cluster QMC re- 
sult. This diagrammatic substitution between the differ- 
ent length-scales in the MSMB method is done by means 
of an analytic Ansatz. 
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One possible implementation of the MSMB method 
which yields a self-energy containing correlations on both 
long and short length-scales is given by the real-space 
Ansatz: 

{^QAfC ~ ^j)' 1^^ ~ ^il - ~2~ (10) 
0, otherwise 

In this formalism, the MSMB self-energy is constructed 
by taking all contributions of lengths up to the lin- 
ear sizes of the small cluster (Nc^'' /2) from the exact 
QMC result and the remaining longer ranged contri- 
butions are complemented by the A-approximated large 
cluster self-energy. Note, the multi-scale self-energy ob- 
tained in Eq. ^| has lost its spatial continuity. Fourier- 
transforming the individual length-scale contributions in 
Eg. [TUI yields 

+ ^ 2 cos{i K2)'t!qIji]{x = i,iujn) 

+ X! 2 cos(j iir2)S^^'^ \x = j,iujn) 
J=-i-+i^-|- 

(11) 

where only diagonal parts of the Fourier transformed self- 
energy are considered. 

J. Hague et. al. in rcf. 3 considered the following, alter- 
native momentum-space Ansatz to combine the different 
length-scales: 

E(^i'')(i^i,zc.„)-s[j^(i^i,zo.„) 

- ^'')(i^i,zc.„) + ^'')(ifi,zc.„) 

(12) 

E(^^^')(i^2,^^„) = Ef^^')(X2,^..„) 

+ ^QMC {K2,lUJn) - E^ [K2,lUJn) 

(13) 

where is the self-energy obtained in the A- 

approximation when implemented on the small cluster. 
The self-energies E'^-'^^ and E^^^ on the small and 
large cluster, respectively, exist on different grid sizes, 
and it becomes necessary to convert self-energies from 
one to the other. This conversion is denoted by a bar 
over the self-energy which denotes an interpolation when 



going from a coarser grid to a finer one and a coarse- 
graining step otherwise. For example, the large cluster 
self-energy E^^^ is constructed from the explicit A- 
approximated self-energy on the large cluster, and two 
interpolated small cluster (denoted by the superscript 

(A^c^^)) self-energies Eg^^^ and S^^"" It is important 
to note that the coarse graining in going from large- to 
small-cluster self-energies is not an averaging over elec- 
tron propagators but the term is used in this context to 
imply an averaging of the self-energy within a cluster-cell. 

While the real space implementation of the Ansatz was 
a straight forward combination of different length-scale 
elements, the momentum implementation interpretation 
is more involved. The A-method provides an estimate for 
the set of self-energy diagrams which convey the long- 
ranged correlations of the system. However, since in re- 
ciprocal space there isn't an explicit separation of length- 
scales, the remaining short length-scale diagrams which 
are to be supplied by the QMC calculation have to be 
identified. This can be accomplished by removing the 
sub-set of A-approximated self-energy diagrams on the 
small cluster from the complete set of the QMC calcula- 
tion hence avoiding a double counting of the correspond- 
ing self-energy contributions. 

Within the traditional form of the DCA, all self- 
energies are inherently causal for each individual cluster. 
One consequence of causality is that — i/mE(A:,w) > 
0. In this Ansatz based implementation of the MSMB 
method however, causality is not inherently guaranteed. 
The combination of self-energy contributions of various 
length-scales in the Ansatz is only assured to yield a 
causal multi-scale self-energy in the limit of the cluster 
sizes approaching one another. Therefore, causality in 
these schemes cannot be guaranteed and hence has to be 
monitored closely throughout. For a more detailed dis- 
cussion of this and the entire momentum-space Ansatz 
see Ref . 0. 

The remaining self-consistent implementation of the 
multi scale method (i.e. two cluster DCA) is similar to 
that of the traditional, single cluster DCA. Fig.[3depicts 
a flow chart of the implementation of both Ansatze in the 
scheme of the overall self-consistency loop. The bottom 
loop shows the already discussed DCA self-consistency 
loop using the QMC as a small cluster solver. For the 
large cluster solver (A-method) the self-consistency is 
similar. In the overall scheme of the Ansatz, the two clus- 
ter problems are combined in a fully self-consistent ap- 
proach: After each iteration of the QMC/DCA loop, the 
corresponding A/DCA contribution on the large cluster is 
evaluated. The Ansatzia used after each step to calculate 
new estimates for the self-energies on both the small and 
large cluster, thus yielding a fully self-consistent solution. 
The corresponding paths are shown in Fig. [S] The con- 
verged Ansatz self-energy will be dominated by the small 
cluster QMC self-energy which contributes the strongest, 
short-ranged correlations. The remaining weaker long- 
ranged correlations are incorporated in the difference of 
A-approximated self-energies between the two clusters. 
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FIG. 5: Flow chart for the Ansatz self-consistent implemen- 
tation of the MSMB/DCA. 



Furthermore, while the coarse-graining (sum over inter- 
nal moment k") is to the small cluster, the self-energy 
in G originates from the large cluster. The resulting 
reducible vertex function F however, is still only de- 
fined for small cluster momenta but incorporates large 
cluster corrections from the coarse grained susceptibility 
X°. Due to these long-ranged contributions, the reducible 
vertex constructed in this manner is not equivalent to a 
QMC evaluated small cluster F. 

The resulting Bethe-Salpeter equation H14() is most eas- 
ily solved for F^ in the following matrix form 



In moving away from a fully self-consistent approach, the 
self-consistency restrictions may be lessened to various 
degrees. Some possible implementations will be discussed 
in further detail in section Hvl 



III. APPROXIMATIONS TO THE VERTEX 

In this section we provide the remaining details of the 
A-approximation. Previous studies, as well as our own re- 
sults, have shown that in the positive U Hubbard model 
the long-ranged contributions to the self-energy are dom- 
inated by the spin and charge fluctuations of the system. 
In contrast, pairing fluctuations are less significant and 
do not change the qualitative multi-scale results, unless 
very low temperatures are considered. The self-energy 
determination in the results section for the lowest order 
in U approximation to considers all three channels. 
However, the remaining higher order cases are restricted 
to only the self-energy contribution of the particle-hole 
channel. An equivalent derivation to the following can 
be trivially extended to the particle-particle channel. 

Starting with the A-approximation introduced in sec- 
tions (Eqs. |H1 and El, the Bethe-Salpeter equation 
(Eq. [7|) for the approximated large cluster reducible ver- 
tex function can be coarse-grained to the small cluster 
(see section 1^ for further details): 

F^{Ki,K[,Q2;iuJn,icj'n,iv) = T^{Ki,K[;iujn,iuj'^) 
xF\KlK[,Q2;iLu';,iu;'„,iiy) 

(14) 

where we coarse grained the bare particle-hole suscep- 
tibility bubble (internal legs in Fig. 0] which are lattice 
Green's functions). We call this coarse grained suscepti- 
bility x° 



(16) 



where the matrix indices correspond to the internal mo- 
menta and frequency. 

At this point, the self-energies can now finally be eval- 
uated using the Dyson equation as illustrated in Fig. O 
On the large cluster this yields: 



T7T)T72) E G{K2 + Q2,i^n + iy) 



x\K[,Q2;iLo',,,iv) X (F"" {K[,M{K2),Q2\^Lo'^,^u;^,^y) 



-A+- 



-ATi 



{K[,M{K2)., Q2\ lio'n, «Wri, iv) - U) 



(17) 



Here we interpolate the small cluster momentum 
M{K2) — > K2 (for details see section IVl|) and subtracted 
U in the parenthesis to prevent an over counting of the 
second order term. 

In the calculation of the transverse spin fluctuation 
part, recall that 2x^ = x^^ where we define x^ ^ the 
correlation function of ct+ and a~ , and x^^ the corre- 
lation function formed from . Then as = 2(x^^ — 



X 



we have that x 



X 



TT 



X 



and F± =fTT-MT. 



This means that for the self-energy on the large cluster. 



T.f'-^'"\K2,l^n) 



A^(l)Ar(2) 



J2 G{K2 + Q2 



(2F^^\i^;, M(if2), Q2; iC^n, + U) X"{K[,Q2; l^'n. 

(18) 



X°{K'l,Q2]lU}n,lv) = - 



N 



k" ,iujn 



X G{K'{ + k" + Q2,iuj[[ + iv) 



(15) 



While for the real-space implementation of the Ansatz, 
knowledge of the large cluster self-energy in the A- 
approximation is sufficient, the momentum-space version 
requires the corresponding self-energy diagrams on the 
small cluster as well. An equivalent calculation on the 
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small cluster yields for the self-energy in A-approximation 



Nc ' K[,Qi,iui'^,u' 

(19) 

where all single-particle propagators have been replaced 
by Gc, including the ones entering the bare bubble x°- 

We want to reiterate that the A-approximated self en- 
ergy thus obtained, is only accurate for long-ranged cor- 
relations. In order to obtain a viable multi-scale solu- 
tion, the neglected short-ranged correlations have to be 
accounted for by means of the Ansatz^ as outlined in the 
previous section. 

A. First Order in U 
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FIG. 6: Second order diagrams for the vertex functions F^^ 
a.) and Fi^ - b.) for an external momentum transfer g = 0. 



The introduced non-perturbative MSMB method re- 
quires a detailed knowledge of the small (QMC) cluster 
irreducible vertex. This evaluation of F however, poses 
a difficult and numerically involved problem. We shall 
therefore initially only consider perturbative approxima- 
tions to the vertex within the MSMB method. In approx- 
imating the irreducible vertex by the frequency indepen- 
dent first order contribution i.e. T{K, K'^iujmiuj'^) = U, 
the self-energy diagrams of Fig. O reduce to the well 
known FLEXiiLiLiS diagrams. In the particle-hole chan- 
nel, the corresponding self-energy is given by 



(20) 



q rn 

The FLEX potential is defined by 



to illustrate some of the properties and problems asso- 
ciated with the Ansatz in momentum space (see Eqs. 1121 
and I13() . After exploring some of the encountered limita- 
tions, we will introduce variations of the original imple- 
mentation which expand the scope of FLEX applicability 
within the MSMB scheme. 



B. Second Order in U 

At lower temperature or for larger U, higher-order 
terms in the vertex are important. The second or- 
der corrections to the irreducible vertex function in the 
A-approximation (i.e. zero external momentum transfer 
(7 = and frequency f„ = 0) are shown in Fig. |B1 The 
irreducible vertex is given in the spin channel by 



V -xAl^^-) 2 1 + xl(g,a.„) 



3 xlhiq.^n) 



(21) 



2 1-X°ft('7,'^«) 



where 



- w„) 



(22) 

Using the FLEX to address the long length-scale prob- 
lem within a multi-scale method was the scope of the 
work by J. Hagu©^. Initially we will return to this simple 
cluster solver which is known to provide qualitatively cor- 
rect long length-scale properties. We will use the FLEX 



■ K'' 



X G{iuJn + iuJni — iuJn" , Ki + K[ — Kl) 

(23) 



and in the charge channel 



X {G{iuJn + il^n' — iLUn",Ki + K[ — K'() 

+ 2 G{iuJn — i^^n' + ioJu-'^Ki - K[ + K'()) . 

(24) 
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C. Full QMC vertex 

In the non-perturbative MSMB approach the full irre- 
ducible vertex, as obtained by the QMC, is considered. 
Within the QMC algorithm we are unable to determine 
r directly and instead calculate the two-particle correla- 
tion function x- The irreducible vertex function can be 
found by inverting the Bethe-Salpeter equation. In the 
spin-channel we have the two particle correlation func- 
tion 

x^ = x" + x°rv (25) 

which is a matrix with elements both in frequency and 
momentum (the same holds for the charge channel). The 
irreducible spin vertex is denoted by and x° is the bare 
spin susceptibility respectively. The correlation functions 
are evaluated in the QMC by sampling over the config- 
uration space (Hirsch-Hubbard-Stratonovich fields^) in 
one of two ways, each one posing its own challenges. 

One possibility is directly evaluating x{i^rn k; itOn', k') 
in frequency space. This requires the individual QMC 
Green's function to be Fourier transformed (FT) from 
the time domain. However, since the calculation is lim- 
ited to finite time intervals At the FT will incur sub- 
stantial high-frequency artifacts. In theory this can be 
improved by means of a high-frequency conditioning, but 
no analytic form is available within the QMC calculation 
to facilitate such a conditioning. This resulting substan- 
tial artifacts are propagated into the irreducible vertex 
function. 

Alternatively x can be evaluated in the time domain 
and only Fourier transformed once the QMC averag- 
ing is complete. This results in an accurate measure- 
ment of the two-particle correlation function by the QMC 
but is significantly more computationally expensive. Al- 
though this approach provides an accurate measure for 
x('''i, 72; Ta, T4), the Bethe-Salpeter equation can not be 
be used to solve for F in the time domain and hence 
the QMC averaged x has to be Fourier transformed first. 
This once again results in similar problems associated 
with the FT. In the high temperature regime where 
higher order corrections in U are not important a per- 
turbative approach to high-frequency conditioning of the 
FT is very successful. However, as the temperature is 
lowered and the significance of higher order corrections 
within the vertex grows, this conditioning results in even 
larger artifacts. Therefore we use the frequency domain 
in our determination of x- 



IV. TECHNICAL ASPECTS 

A MSMB approach based on a first order approximated 
vertex is similar to the limited FLEX-hybrid approach 
previously considered by J. Hague et. aim-. In their work 
the intermediate length-regime was addressed using the 
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FIG. 7: Imaginary part of the self-energy at lowest Mat- 
subara frequency as obtained by various cluster solvers and 
FLEX/MSMB method at /3 = 8, C/ = = 1.0, and n = 0.75. 
The MSMB cluster sizes are ivi^' = 8 and A^i^^ = 32. 



FLEX which was incorporated within a multi-scale ap- 
proach using the fully self-consistent momentum space 
Ansatz (sec Eas. 1121 and However, in the regime of 
stronger couplings and/or lower temperatures, where a 
significant contribution of longer ranged correlations is 
to be expected, the method fails. In this section we pro- 
pose changes to the implementation of the multi-scales 
method which significantly improves its range of applica- 
bility by 1) restriction of the Ansatz to the large cluster, 
2) modification of the cluster conversion, and 3) removal 
of the self-consist implementation of the Ansatz on the 
large cluster. 

In the original inception'^ , the Ansatz is evaluated fully 
self-consistently within both the QMC and FLEX as dis- 
cussed in section III Dl The self-energy contributions of 
these two clusters solvers are linked by the momentum 
space Ansatz. However, long ranged-correlations, as are 
described by the large cluster FLEX, are assumed weak. 
Hence, their presence in the effective medium of the small 
cluster only minutely effects the QMC self-energy. We 
thus neglect the Ansatz on the small cluster (Eq. I12|l . 
With this modification, the small cluster problem can be 
solved independently of the large cluster problem. This 
results in a significant reduction in numerical complexity. 

With the removal of the small cluster Ansatz condi- 
tion, only a unidirectionally cluster conversion of small 
cluster self-energies to the large cluster remains. In this 
work, a periodic cubic spline interpolation^^ is employed. 
This provides a good approximation for the multi-scale 
self-energy in the high temperature/small U limit. Fig.[7| 
shows the imaginary part of the exact, QMC calculated 
and FLEX self-energies on both clusters calculated at 
high T. In this regime, correlations are short ranged in 
nature and thus well described by the small cluster itself, 
while remaining long-ranged features, which are provided 
by the large cluster FLEX, are insignificant. The interpo- 
lated small cluster self-energy accurately replicates that 
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FIG. 8: Flow chart depicting two independent self-consistent 
DCA calculations combined via an Ansatz to construct a 
MSMB self-energy. 



of the explicit large single cluster QMC result. This re- 
sults in a vanishing FLEX contribution within the Ansatz 
as the large and interpolated small cluster FLEX results 
are identical and thus cancel each other (see superim- 
posed, lower set of curves in Fig.|7|). 

In this high temperature regime (e.g. Fig. [7|l the fully 
self-consistent Ansatz as well as our modified one re- 
main numerically stable. However, one inherent lim- 
itation of the approach is already apparent. The ob- 
served magnitude of the negative imaginary part of the 
QMC self-energy is significantly smaller than the FLEX 
result. This overestimation of the self-energy by the 
FLEX deems a first order approximation to the vertex 
to be insufficient. A self-consistent implementation of 
the Ansatz aggregates the problem as the Ansatz self- 
energy used to initialize each FLEX iteration is similar 
in magnitude to the QMC self-energy. This results in an 
insufficient damping of the FLEX potentials and the sub- 
sequent overestimation of the FLEX self-energy renders 
the large cluster self-energy calculation numerically un- 
stable. This is an inherent problem in the self-consistent 
approach. In an attempt to enhance the numeric stabil- 
ity of the MSMB method we remove the self-consistent 
implementation of the Ansatz on the large cluster. This 
leaves two independent DCA calculations for the small 
and large cluster which, after the individual problems are 
converged, are combined by the Ansatz. The resulting 
self-consistency scheme is depicted in Fig. |H1 However, 
the effective media embedding the two cluster problems 
now are unaware of correlations as determined by the 
other cluster solver technique. Therefore, the effective 
medium of the large cluster now lacks the explicit short 
ranged correlations as provided by the QMC. Similarly, 
the effective medium of the small cluster only contains 
long-ranged correlations as provided by the dynamical 
mean-field approximation. 

The third modification enters in the determination of 

(^(1)-) 

the small cluster FLEX self-energy ^pf^x- Rather than 
calculating it explicitly, we obtain it by coarse graining 
the large cluster self-energy onto the small cluster. The 
combination of all changes introduced up to this point 
significantly increase the range of low temperatures which 
can be addressed by a perturbative MSMB approach. 

Throughout the remainder of the paper we will restrict 
ourselves to the non-self-consistent momentum space 
Ansatz as outlined in this section. It indeed succeeds 
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FIG. 9: Imaginary part of the self-energy at lowest Matsub- 
ara frequency as obtained by the MSMB method using first 
(FLEX) and second order approximated irreducible vertices 
F in comparison to the large single-cluster QMC results at 
/3 — 31, U = W = 1.0, and n = 0.75. Multi-scale results are 



for cluster sizes iVi^' = 8 and Ni; 



r(2) 



32. 



in addressing some of the problems encountered in the 
original implementation and yields a numerically stable 
MSMB method. 



V. RESULTS 

Throughout the remainder of the paper, we evaluate 
the quality of the MSMB method by comparing the self- 
energy to that of a 32 site, single cluster DCA/QMC 
calculation. 



A. First Order in U 

We begin the evaluation of the MSMB method by 
considering the first order approximation to the irre- 
ducible vertex. The previous section introduced a nu- 
merically stable approach to the perturbative MSMB 
method. This allows for the FLEX based MSMB treat- 
ment of stronger coupling/lower temperature regimes. 

The major limitation of this method however remains 
in the still significant difference in magnitude of the QMC 
and FLEX self-energies leading to an overestimation of 
long length-scale features introduced by the FLEX. In 
Fig-Elit is quite apparent that at lower temperatures the 
FLEX MSMB implementation overestimates the size of 
the long length-scale features i.e. amplitude of the oscil- 
lations in the imaginary part of the self-energy. This is 
yet further indication that a bare approximation to the 
vertex is inadequate to address the intermediate length 
regime. 
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B. Second order in U 

We move to include second order corrections in the 
vertex which is expected to mitigating the effects of the 
underestimation to the vertex in the FLEX. This lessens 
the difference in magnitude between the QMC and A self- 
energy within the Ansatz (as observed in Fig. 01 and thus 
increase the applicability of the MSMB method to lower 
temperatures. 

Fig. 1^1 compares the imaginary part of the self-energy 
at lowest Matsubara frequency for various degrees of ap- 
proximation. In contrast to the gross overestimation of 
the first order approximation to F, the inclusion of sec- 
ond order in U corrections within the MSMB method 
successfully captures the long length-scale features of the 
large single-cluster QMC self-energy throughout most of 
the Brillouin zone. The largest deviations in the multi- 
scale self-energy are found about the corners of the Bril- 
louin zone i.e. k = ±7r. Within the Ansatz, the differ- 
ence of the small and large cluster A self-energy provides 
the long length-scale features thus partially restoring the 
self-energy information which was lost by coarse grain- 
ing to the small QMC cluster. However, using a A self- 
energy with a negative imaginary part which is consid- 
erably larger than that of the small cluster QMC result 
in an overestimation of this correction, predominately in 
regions encompassed by steep gradients in the self-energy. 

Fig- El further illustrates the pathology of the Ansatz 
associated with the difference in magnitude of the clus- 
ter self- energies. We show the imaginary part of the self- 
energy as obtained by the various A-approximations on 
the large cluster in comparison to the large single-cluster 
QMC result. The magnitude of the negative imaginary 
part of the self-consistent FLEX (first order in U approx- 
imation) self-energy is, as was previously indicated, very 
large compared to the magnitude of the QMC result. The 
new approach of including second order corrections in U 
in the vertex function, succeeds in yielding a self-energy 
which resembles that of the QMC more closely. Similarly 
to the numeric instability encountered with the FLEX 
approximation, the self-consistent determination of the 
A self-energy using the large cluster Ansatz also yields a 
larger self-energy but doesn't encounter the catastrophic 
divergence of the FLEX. We find that a second order F 
combined with a non-self-consistent Ansatz provides the 
best multi-scale solution to the problem. 

Now that we have established the viability of the 
MSMB technique we take a closer look at the depen- 
dence of the small cluster size on the quality of the multi- 
scale result. Within the MSMB method the QMC consti- 
tutes the computationally most expensive part. Hence, 
we want to restrict the calculations to the smallest pos- 
sible QMC cluster {Ni^^) without any significant loss of 
quality in the multi-scale results. It is therefore impor- 
tant to study the dependence of the multi-scale results 
on the size of the small cluster. 

Fig. ^] shows the multi-scale self-energy for a variety 
of small cluster sizes as obtained by the MSMB method 
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FIG. 10: Imaginary part of the A-approximated self-energy at 
lowest Matsubara frequency for the large cluster (A^i'^^ — 32) 
using various cluster solvers in comparison to the single cluster 
QMC result at /3 = 31, f/ = W = 1.0, and n = 0.75. Also 
shown is the overestimated Ansatz self-consistent self-energy 
of the MSMB method for the second order approximated P. 



using the second order approximation to the A vertex. It 

is apparent that a cluster size of Np'^ = 4 is too small 
to adequately capture the main k dependence (overall 
magnitude) of the QMC's self-energy. Considering such 
a small QMC cluster grossly misrepresents the range of 
the QMC self-energy and is not recovered in the MSMB 
method. This is an indication that correlations beyond 
the length-scales of the small cluster are still significant 
and not sufficiently well approximated by the A-method. 
A small cluster size of ivi"'^' = 8 on the other hand, ap- 
pears adequate and only slightly underestimates the self- 
energy in the vicinity of A: = ±7r.In this area we continu- 
ally observe significant remnants of the coarse graining of 
the QMC which is inadequately restored by the MSMB 
method. 

Up to this point we have focused our investigation on 
the momentum dependence of the self energy at small 
Matsubara frequency. Since the visual representation of 
the self-energy features at multiple, larger Matsubara fre- 
quencies would be cumbersome, we proceed to further 
illustrate the strengths of this MSMB technique by fo- 
cusing on the presence of spin-charge separation in the 
system which is manifest in the full frequency dependent 
self-energy. 

One dimensional systems have been shown to be non- 
Fermi liquids. Amongst other unique features they are 
known to exhibit spin-charge separationi2Si2iiS2i. This 
very intriguing property manifests itself by the com- 
plete decoupling of spin and charge degrees of freedom. 
The single-particle spectra of such systems exhibits two 
unique peaks corresponding to either spin or charge exci- 
tations which move independently of each other. Due to 
the involved nature of extracting spectra from the MSMB 
method, we are at this time unable to directly identify 
the presence of any such feature in our results. One char- 
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FIG. 11: Imaginary part of the self-energy at lowest Matsub- 
ara frequency as obtained by the MSMB approach using the 
second order A-approximation and large single-cluster QMC. 
Results are for various small cluster sizes iVi^' at /3 = 31, 
V = W = 1.0, and n = 0.75. 



acteristic of such a separation however, is the presence of 
two distinct velocities corresponding to charge and spin 
respectively. Following an approach by Zacher^^, we ex- 
amine the MSMB Matsubara frequency Green's function 
to the possible presence of spin-charge separation. This 
is done by fitting Green's function obtained from the 
MSMB method to that of the Luttinger model solution 



q{.LM) , X 



(26) 



where v\ and are the spin and charge velocities and c 
is a normalization constant. We use the approximation 
for the correlation exponent Kp = 1 which is deemed suf- 
ficient by Zacher for purposes of identifying the presence 
of two different velocities. In order to accurately perform 
a fit with our data we need to additionally coarse grain 
(see Eq. O the Luttinger liquid Green's function to ob- 
tain a fitting function in-line with the idea of the DCA. 
This fit yields values for both the spin and charge veloci- 
ties when fitted at fc 7r/2 which is the DCA momentum 
closest to the Fermi wave vector. This specific choice of 
k value is motivated by fact that the Luttinger model 
solution is based on a low energy approximation of the 
Hubbard model where a linearized dispersion around the 
Fermi vector is assumed. For the parameters in Fig. 1121 
the Fermi wave vector (kp « 1.2) falls into the cluster 
cell about fc = 7r/2. 

Fig. 1 121 shows the spin and charge velocities {v\ and V2 
respectively) obtained by using the single cluster QMC 
results for smaller cluster sizes and multi-scale results for 
larger ones. For intermediate cluster sizes, where both 
the QMC and the MSMB method are feasible, we see a 
good match between the two methods. For larger cluster 
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FIG. 12: Spin and charge velocities [v\ and v2 respectively) 
obtained by fitting the different results with the Luttinger 
Green's functions (Ea. l26ll about cluster momentum fe = 7r/2 
for /3 = 31, [/ = = 1.0, and n = 0.75. Multi-scale results 
are for a small cluster size hS'^p = 8 and different large cluster 
sizes, and the QMC velocities were obtained from a single 32 
site cluster calculation. 



sizes where the QMC approach becomes infeasible the 
multi-scale results fall close to the QMC extrapolated 
values. Furthermore, as we extrapolate the multi-scale 
results to the infinite cluster size limit we find two dif- 
ferent velocities: v\ — 0.311 and V2 = 0.674. These 
values compare rather well to those obtained by Zacher 's 
et. alm^ grand-canonical QMC calculation for a cluster 
size of 64, /? = 80 and all other parameters equal to ours: 
Vl ^ 0.293 ± 0.019 and V2 = 0.513 ± 0.023. Although 
the temperature in Fig. 1121 is slightly higher compared to 
that of Zacher's results, we observed only a small tem- 
perature dependence of the fitted velocities in our calcu- 
lations. The similarity between the two results arc quite 
remarkable considering that our calculations were based 
on a substantially smaller 8 site QMC cluster and hence 
less subject to the sign-problem . These are yet further 
indications that the MSMB method indeed successfully 
captures the long length-scale physics of the model and is 
in good quantitative agreement with large single-cluster 
QMC calculations. 



C. Full QMC Vertex 

The previous two sections considered initially a first or- 
der approximated irreducible vertex function (i.e. FLEX) 
and then proceeded to show the significant improvement 
one can achieve by including second order corrections to 
r within the MSMB method. Both of these perturba- 
tive approaches however resulted in an overestimation of 
the long length-scale features of the self-energy at low 
temperatures. We now return to the non-perturbative A- 
approximation using the irreducible QMC vertex. This 
approximation is expected to significantly improve the 
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FIG. 13: Imaginary part of the self-energy at lowest Mat- 
subara frequency as obtained by the MSMB method using 
second order approximated and full QMC evaluated Fs in 
comparison to a large single cluster QMC results at /3 = 31, 
U — = 1.5, and n — 0.75. Multi-scale results are for 
cluster sizes iVi^' = 8 and = 32. 



results, especially at lower temperatures where higher 
order corrections to the bare vertex are important and 
measurably contribute to the self-energy. 

Fig . 1 1 31 compares the imaginary part of the multi-scale 
self-energies obtained by considering both the second or- 
der vertex and the full QMC vertex, to the numerically 
exact large single-cluster QMC self-energy. While the 
resulting multi-scale self-energy utilizing the full vertex 
function gives a better approximation for some momenta, 
the second order vertex solution remains preferable for 
others. This behavior may be somewhat unexpected 
since the inclusion of higher order corrections in the ver- 
tex was expected to further improve the overall quality 
of the multi-scale solution. However, we believe that this 
discrepancy arises from difhculties in extracting the ex- 
act vertex related to problems with high frequency con- 
ditioning in the QMC and is not an intrinsic problem 
of the method. As a consequence, the consideration of 
second order corrections in the A-approximation remains 
the most useful at this time. It successfully allows for the 
exploration of lower temperatures - a regime where the 
reducible vertex function develops a richness in features. 
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FIG. 14: Imaginary part of the self-energy at lowest Mat- 
subara frequency as obtained by the MSMB method using 
the second order approximated F in conjunction with the 
real-space and momentum-space based Ansatz at /3 = 31, 
U = W — 1.0, and n = 0.75. Also shown is the large single- 
cluster QMC result. Multi-scale results are for cluster sizes 
iVi^' = 8 and iVi^' = 32. 



and will hence be used to investigate this Ansatz imple- 
mentation as well. It should be noted that both versions 
of the Ansatz are treated identically as far as the imple- 
mentation of self-consistence is concerned (see Sec. lV A|l . 

Fig. El compares the multi-scale self-energy of the two 
Ansatz implementations to the large single cluster QMC 
result. It is quite apparent that the real-space implemen- 
tation of the MSMB method is inferior. The momentum- 
space Ansatz provides a self-energy more closely resem- 
bling the exact result throughout the Brillouin zone. 

The real-space Ansatz provides an intuitive and simple 
way to combine the different length-scales of the problem 
in contrast to a more complicated implementation in mo- 
mentum space. Although inferior, the real-space Ansatz 
remains a viable, numerically stable alternative. It fur- 
thermore provides an alternative means of interpolating 
the small cluster self energy by neglecting the large clus- 
ter self-energy contributions in Eq. 1111 which were found 
to be negligible in this approach. The resulting interpo- 
lated small cluster QMC self-energy (not shown) closely 
resembles that of the the real-space Ansatz based MSMB 
method. 



D. Real Space Ansatz 

Up to this point, the various approximations to the ir- 
reducible vertex function within the MSMB technique 
were implemented using the momentum-space Ansatz. 
We now return to the previously introduced implemen- 
tation of the real-space Ansatz (see Eg. Ill|) . The signif- 
icant difference is that the separation of length-scales in 
the real-space approach is straightforward, and therefore 
over-counting of diagrams is not an issue. The second 
order vertex approach was shown to be widely successful 



VI. NUMERICAL CONSIDERATIONS AND 
OUTLOOK 

In the development of this MSMB method we were 
forced to employ various approximations due to current 
computational limitations. The largest numerical con- 
cession was made in section Hill where we restricted the 
calculations of F'^ to the small cluster (see Eq. I14|l . Ide- 
ally, the irreducible vertices F would be interpolated and 
the full reducible vertex evaluated on the large cluster in 
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FIG. 15: Parquet equation relating the transverse particle- 
hole irreducible vertex F^^^j to the fully irreducible vertex A 
plus contributions from the longitudinal and particle-particle 
cross channels. Similar relations apply for the remaining 
channels (not shown). 
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FIG. 16; Lowest order non-local corrections to the fully- 
irreducible vertex A (left) and the vertex F (right). 



turn. This calculation however would scale as (TVi Ni)'^ 
{Ni is the number of time-slices in the QMG^) and hence 
provide little advantage over a single cluster QMC cal- 
culation. We therefore restrict the evaluation oi F^x'^ in 
Eq. ^] to the small cluster and interpolate the product 
of and F'*' to the large cluster. 

With the onset of peta-scale computing we will be able 
to make two fundamental improvements to the MSMB 
approach in the near future. Initially, we will gain the 
ability to include the fully momentum and frequency de- 
pendent F in our calculation, thus eliminating the ne- 
cessity of the A-approximation (Eq. 15)). Inherent in this 
modification is an explicit account of the correct short 
ranged physics hence removing the need of the Ansatz. 
However, the memory and CPU requirements for this 
type of calculation scale as (Nc^^^ Ni)^. For a large clus- 
ter with Nc = 16 and Ni = 100, this would require 66G 
of double precision complex storage, far exceeding the 
memory associated with a single CPU. These stagger- 
ing memory requirements can currently only be met by 
some shared-memory parallel processing (SMP) super- 
computers. 

In the second improvement, the approximation of the 
large cluster F by the small cluster QMC one can be 
replaced by one utilizing the fully irreducible vertex A 
(the vertex which is two-particle irreducible in both the 
horizontal and vertical plane). This results in a self- 
consistent renormalization of F via the Parquet equa- 
tiongi2*S4, and hence an inclusion of long ranged corre- 
lations in the crossing channel (see Fig. I15|l which are 
missing in both the A-approximation and the F-based 
approximation described above. 

The superiority of the latter approach becomes clear in 
the high-dimensional limit, where it is A, not F, which 
becomes local. This can be shown by considering the 
simplest non-local corrections to the respective vertices 
A and F in Fig. El The boxes represent a set of graphs 
restricted to site i (local) and j (neighboring) respec- 
tively. In the limit of high dimensions, each site i has 2D 
adjacent sites j. The contributions of each leg within the 
vertex in the limit Z) — )■ cxd is G(r) ~ (for details 

see Ref. @) . This results in a contribution to the correc- 
tion of O(D-i) for the two legs in F and 0{D-^/^) for A. 



Thus, the non-local corrections to A including all neigh- 
boring sites j falls off as D~^/'^ and becomes local in the 
infinite-dimensional limit. In contrast, the corrections to 
F remain of order one. Therefore, in the high-dimensional 
limit, A is local while F has non-local corrections. In fi- 
nite dimensions, we would expect that A is more compact 
than F whenever the single-particle Green function falls 
quickly with distance. Then A should always be better 
approximated by a small cluster calculation than F. (De- 
spite the fact that F has non-local corrections, one can 
easily show that in the high-dimensional limit, all of the 
methods discussed here will yield the same self energy 
and susceptibilities since the non-local corrections to F 
fall on a set of zero measure points). 

In employing the solution to the Parquet equation in 
a MSMB method we would be able to resolve two ma- 
jor limitations of the current approach: 1) An imple- 
mentation considering the full frequency and momentum 
dependent vertex will be devoid of the causality prob- 
lems associated with the self-energy mixing of the two 
cluster sizes. 2) The approach constitutes a conserving 
approximation for the large cluster self-energy. Given 
these potential gains of a future method, we have to 
stress the extensive computational demands associated 
with this approach. While in a F based implementation 
a trivial numerical parallelization of the problem leaves 
manageable demands, the complex nature of the Parquet 
approach requires substantial future development. 



VII. CONCLUSION 

We have introduced a numerically feasible MSMB ex- 
tension to the DCA. In this method the lattice problem is 
mapped onto that of two embedded clusters, dividing the 
problem into three length-scales. Correlations on each of 
the length-scales are approximated commensurate with 
the strength of the correlation on the respective scale. 
The intermediate length regime, which bridges the ex- 
plicit treatment of short ranged correlations by means of 
the QMC to the long ranged dynamical mean-field one, 
is addressed in a diagrammatic long-wave length approx- 
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imation based on the two particle irreducible vertex of 
the small cluster. The first order approximation to the 
vertex results in the FLEX but including higher order 
corrections result in a substantial better multi-scale re- 
sult when compared to explicit large cluster QMC cal- 
culations. This can be attributed to the significance of 
higher order diagrams at lower temperatures. We pro- 
ceeded to show that our MSMB results indicate spin and 
charge separation and the obtained velocities compare 
favorably to significantly larger finite size QMC calcula- 
tion. The inclusion of the explicit QMC calculated vertex 
is currently still limited but further work in this direction 
looks promising. However, in any of the introduced im- 
plementations, the MSMB approach provides a means to 
adequately address large cluster problems on all length- 
scales at significant lower computational expense. 
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